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ON  THE  THEORY  OF  OSCILLATING  ALRFOUS  OF  FINITE  SPAN 
IN  SUBSONIC  COMPRESSIBLE  FLOW 
By  Eric  Eeissner 


SUMMARY 


The  problem  of  the  oscillating  lifting  surface  of  finite  span  in 
subsonic  compressible  flow  is  reduced  to  an  integral  equation.  The 
kernel  of  the  integral  equation  is  approximated  by  a  simpler  expres¬ 
sion,  on  the  basis  of  the  assumption  of  sufficiently  large  aspect 
ratio.  With  this  approximation  the  double  integral  occurring  in  the 
formulation  of  the  problem  is  reduced  to  two  single  integrals,  one  of 
which  is  taken  over  the  chord  and  the  other  over  the  span  of  the 
lifting  surface.  On  the  basis  of  this  reduction  the  three-dimensional 
problem  appears  separated  into  two  two-dimensional  problems,  one  of 
them  being  effectively  the  problem  of  two-dimensional  flow  and  the 
other  being  the  problem  of  spanwise— circulation  distribution.  Earlier 
results  concerning  the  oscillating  lifting  surface  of  finite  span  in 
incompressible  flow  are  contained  in  the  present  more  general  results. 


INTRODUCTION 


The  present  report  is  concerned  with  the  problem  of  the  oscillating 
airfoil  of  finite  span,  within  the  frame  of  the  linearized  lifting- 
surface  theory.  The  aim  of  this  study  is  the  development  of  a  theory 
which  incorporates  simultaneously  the  effects  of  three-dimensionality 
of  the  flow  and  of  compressibility  of  the  fluid.  As  an  exact  solution 
of  this  problem,  even  within  the  limitations  of  the  linearized  theory, 
presents  very  great  difficulties,  it  is  worth  while  to  work  toward  an 
approximate  theory  which  is  valid  provided  the  aspect  ratio  of  the 
lifting  surface  is  not  too  small. 

The  author  has  previously  obtained  results  of  this  nature  for  the 
case  of  incompressible  flow  (references  1  and  2).  In  this  earlier 
work  the  known  results  for  the  problem  of  two-dimensional  incompressible 
flow  were  contained  as  a  special  case.  The  present  work  generalizes 
these  results  so  as  to  take  account  of  compressibility  in  the  subsonic 
range.  Thus  the  results  of  this  report  consist  of  a  system  of  equations 
which  contain  as  special  cases  both  the  author’s  results  for  the  wing 
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of  finite  span  in  incompressible  flow  and  the  results  of  Possio*  s  theory 
of  two-dimensional  compressible  flow  (reference  3) • 

The  scope  of  the  present  results  may  briefly  be  described  as 
follows.  The  starting  point  of  the  work  is  an  integral-equation 
formulation  of  the  problem  of  the  lifting  surface  of  finite  span. 

The  integrals  which  occur  are  double  integrals  and  the  functions  to 
be  determined  are  functions  of  two  independent  variables.  The  essential 
step  of  the  present  work  is  to  replace  the  actual  kernel  of  the  integral 
equation  by  an  approximate  kernel  in  such  a  way  that  the  double  integrals 
are  reduced  to  single  integrals  over  the  range  of  either  one  of  the  two 
independent  variables .  In  this  way  the  problem  is  'reduced  to  two 
problems  which  are  to  be  solved  separately.  The  first  of  these  two 
problems  is  of  the  same  nature  as  the  Possio  problem  of  two-dimensional 
compressible  flow.  The  second  of  these  problems  is  of  the  same  nature 
as  the  problem  of  the  Prandtl  lifting— line  theory  for  the  wing  of 
finite  span  in  uniform  motion. 

As  in  the  theory  of  incompressible  flow,  this  reduction  of  the 
double— integral  problem  to  two  single— integral  problems  depends 
crucially  on  the  assumption  of  sufficiently  large  aspect  ratio.  While 
“sufficiently  large"  aspect  ratios  might  be  thought  to  be  aspect  ratios 
of  about  three,  definite  statements  of  this  nature  must  be  based  on 
experimental  evidence,  as  long  as  no  exact  solutions  exist  for  the 
three-dimensional  problem  of  the  oscillating  lifting  surface  in  com¬ 
pressible  flow. 

It  is  perhaps  worth  while  to  state  explicitly  that  the  present 
problem  is  quite  different  from  the  corresponding  problem  for  supersonic 
flow. 


It  might  also  be  added  that  there  are  reasons  to  believe  that  it 
is  not  satisfactory,  even  approximately,  to  superimpose  aspect— ratio 
corrections  for  incompressible  flow  and  compressibility  corrections 
for  two-dimensional  flow  in  order  to  obtain  corrections  for  the  combined 
effect.  This  latter  point  is  one  of  the  reasons  for  the  present  study. 

In  this  first  report  the  work  is  carried  to  the  point  where  the 
double— integral  equation  is  reduced  to  an  equation  containing  single 
integrals  only  over  the  quantities  to  be  determined.  Further  develop¬ 
ments  will  be  given  in  a  subsequent  report. 

This  work  was  conducted  at  the  Massachusetts  Institute  of  Technology 
under  the  sponsorship  and  with  the  financial  assistance  of  the  National 
Advisory  Committee  for  Aeronautics. 
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SYMBOLS 

Cartesian  coordinates 

main-stream  velocity  in  X-direction 

time 

defined  by  equation  of  lifting  surface  Z  =  H(X,Y,t) 

components  of  velocity  change  caused  by  presence 
of  lifting  surface 

region  in  X,Y-plane  occupied  by  projection  of 
lifting  surface 

density  of  stream  flowing  with  velocity  U 

density  and  pressure  changes,  respectively, 

associated  with  velocity  changes  u,  v,  and  w 

velocity  of  sound  in  main  stream  (a2  =  dp/dpQ) 

coordinate  of  trailing  edge  of  E0 

Cl 

potential  of  velocity  changes  u,  v,  end  v 
circular  frequency  of  oscillation 
Mach  number  of  main  stream  (M  =  U/a) 
real  part  of 

a  length  to  be  identified  with  the  semichord 
of  Ba  at  midspan 

reduced— frequency  parameter  (k  =  CDb/U) 
dimensionless  coordinates  defined  by  equation  (19) 
function  defined  by  equation  (21) 


parameter  defined  as  p  =  — — — 

1  -  M2 
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parameter  defined  as  k  =  - ^ 

k 

parameter  defined  as  v  -  — 

1  -  M2 

region  in  x,y— plane  corresponding  to  region  Ra 
in  X,Y-plane 

coordinate  of  trailing  edge  of  R&* 

coordinate  of  leading  edge  of  Ra* 

region  in  x,y— plane  consistixg  of  the  strip  to  the 
right  of  the  trailing  edge  of  Ra* 

entire  x,y-plane  except  for  regions  Ra*  and  R** 

function  defined  by  equations  (36) 

function  defined  by  equation  (37) 

variable s  of  integration  in  accordance  with 
equation  (40) 

defined  by  equation  (hi) 

Hankel  functions  of  second  kind,  and  of  zeroth  and 
first  order,  respectively 

auxiliary  variable  of  integration 

function  defined  by  equation  (51) 

function  defined  by  equation  (52) 

functions  defined  by  equations  (5*0  j  (55) >  (56) 

order  of  magnitude  of 

function  defined  by  equation  (38) 

function  defined  by  equations  (73)  and  (83) 
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cr 

auxiliary  variable  of  integration 

T 

auxiliary  variable  of  integration 

* 

a 

function  defined  by  equation  (76) 

S 

ratio  of  semispan  to  semichord  at  midspan 

Tj* 

defined  as  tj*  =  rj/s\/l  —  M2 

b* 

local  semichord  divided  by  semichord  b  at  midspan 

a  quantity  indicating  amount  of  sweep  and  defined  by 
equation  (86) 

X* 

dimensionless  coordinate  defined  as  x*  =  (x  —  Xnp'jfb* 

K*,v*,n* 

quantities  defined  by  equation  (87) 

g*,x* 

quantities  defined  aB  g*(x*,y*)  =  g(x,y);  X*(x*,y*)  =  X(x,y) 

N 

function  defined  by  equation  (89) 

SM 

function  defined  by  equation  (90) 

THE  BOUNDARY— YALUE  PROBLEM  OF  THE 
OSCILLATING  LIFTING  SURFACE 


It  is  assumed  that  a  nearly  plane,  impenetrable  surface  is  put  into 
the  path  of  an  inviscid  flowing  fluid  which,  except  for  the  effect  of 
this  surface,  possesses  a  uniform  velocity  U  in  the  direction  of  the 
positive  X-axis.  The  impenetrable  surface,  henceforth  called  lifting 
surface,  is  taken  to  lie  nearly  in  the  X,Y— plane  and  its  equation  is 
written  in  the  form  Z  =  H(X,Y,t).  When  H  =  0  no  disturbance  is 
caused.  When  the  lifting  surface  is  not  exactly  plane  and  parallel  to 
the  direction  of  U  the  velocity  components  (U,0,0)  are  changed 
into  (U  +  u,v,w)  where  u,  v,  and  w  depend  on  the  form  of  the 
function  H  and  on  the  shape  of  the  region  Ea  which  is  the  projection 
of  the  lifting  surface  onto  the  X,Y— plane. 

The  disturbances  caused  by  the  presence  of  the  lifting  surface 
are  assumed  to  be  small,  in  the  sense  that  the  differential  equations 
and  boundary  conditions  of  the  problem  are  linearized  with  respect  to 
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the  disturbance  velocity  components  u,  v,  and  w  and  with  respect  to 
the  pressure  and  density  changes  p  and  '  p  caused  by  the  presence  of 
the  lifting  surface . 

Under  these  conditions  the  differential  equations  of  the  problem 
are  the  following: 


■  du  +  -g  du  _  d  /P  \ 

at  ax  ”  axlpo/ 


(1) 


U  —  =  -  d  /P  \ 

ax  az\p0y 


p  =  a2p  (5) 

The  quantity  pQ  in  equations  (l)  to  (4)  is  the  density  in  the  fluid 

flowing  without  disturbance  and  the  quantity  a  in  equation  (5)  is  the 
velocity  of  sound  in  the  undisturbed  fluid,  that  is,  a2  =  dp/dpQ. 

The  boundary  condition  of  no  relative  normal  flow  at  the  lifting 
surface  is  satisfied,  within  the  frame  of  the  linearized  theory,  instead 
of  on  the  lifting  surface  itself,  on  the  projection  of  this  surface  onto 
the  X,Y— plane. 


X,Y  inside  Ea,  w  =  +  U  ^  (6) 

The  form  of  condition  (6)  (which  holds  on  both  sides  of  the  lifting 
surface)  indicates  that  w  is  an  even  function  of  Z.  From  equation  (3) 
it  follows  then  that  p  is  an  odd  function  of  Z  and  the  condition 
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that  the  pressure  disturbance  p  is  continuous,  except  when  passing 
across  the  lifting  surface,  means  that  for  Z  =  0  and 

X,Y  outside  Ra,  p  =  0  (7) 

On  the  basis  of  conditions  (7)  and  (6)  the  problem  may  be  considered 
to  consist  in  the  determination  of  u,  v,  w,  and  p  in  the  half 
space  Z  ^  0  with  the  boundary  Z  =  0. 

In  addition  to  conditions  (6)  and  (7)  the  following  further 
conditions  are  prescribed  in  order  to  obtain  an  ■unambiguous  solution. 

At  the  trailing  edge  of  the  lifting  surface, 

X  =  Xp(Y),  p  is  finite  (8) 

Finally,  it  is  postulated  as  a  condition  "at  infinity"  that  energy-  is 
traveling  outward  without  reflection,  in  a  manner  to  be  defined  more 
precisely  in  what  follows  for  the  case  of  simple  harmonic  motion. 


VELOCITY— POTENTIAL  FORMULATION  OF  THE  PROBLEM 


It  can  be  shown  that  the  problem  may  be  solved  as  stated  in 
equations  (l)  to  (8)  by  means  of  a  velocity  potential  0,  in  terms 
of  which 


Combination  of  equations  (9),  (l),  (2),  and  (3)  results  in  the  following 
expression  for  the  pressure  change  p: 
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Combination  of  equations  (10),  (5),  (9),  and  (4)  results  in  the 
following  differential  equation  for  0: 


Equation  (ll)  is  to  he  solved  in  the  half  space  Z  >  0  subject  to  the 
following  boundary  conditions  at  Z  =  0: 

X,Y  inside  Ea,  g  -  g  +  °  g  (l2) 

X,Y  outside  Ra,  ^  +  U  ^  =  0  ( 13 ) 

ot  oX 

X  =  Xt(y)j  ^  +  U  ^  is  finite  (l4) 

1  St  dX 

and  subject  to  the  condition  of  no  energy  reflection  "at  infinity." 

In  what  follows  attention  is  restricted  to  the  case  of  simple 
harmonic  motion  and  the  following  equation  is  written: 


0(X,Y,Z,t)  =  eia)t0(X,Y,Z) 


(15) 


with  corresponding  expressions  for  H  and  p. 


1It  is  perhaps  not  entirely  superfluous  to  indicate  that  this  is 
meant  in  the  sense  that,  corresponding  to  a  surface  equation  Re(eia^H), 
there  is  a  pressure  distribution  Re(eittytp). 
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Equation  (ll)  now  assumes  the  form. 


d20  a2H 

ax2  +  ay2  +  az2 


-  M1* + D  s)2?  - 


(16) 


Equation  (10)  becomes 


P  =  “Po(^0  +  U^j 


(17) 


How  the  following  dimensionless  variables  and  parameters  are 
introduced: 


M  = 


k  = 


cub 


(18) 


X  =  b 


=  \jl-U‘ 


2  Z 
b 


(19) 


The  differential,  equation  (l6)  then  assumes  the  form 
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Equation  (20)  is  reduced  further,  for  the  purpose  of  eliminating 
first-derivative  terms,  "by  the  following  substitution: 


0  =  e 

Choose 


(21) 


(22) 


and  obtain  as  the  equation  for  ty, 

+  K^Ir  =  0  (23) 


where 


K 


m 

1  -  M2 


(2k) 


Now  p  and  the  boundary  conditions  must  be  expressed  in  terms  of 
the  new  independent  variables  x,  y,  and  z  and  in  terms  of  the  new 
dependent  variable  i|r. 

From  equations  (17)  and  (19),  it  follows  first  that 


'5  =  -^  +  l)  (25) 

Combination  of  equations  (25)  and  (21)  gives 

P:_M.^vt  +  D 


(26) 
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where  V  is  defined  as 


The  boundary  condition  (12)  becomes 

x,y  inside  V.  g  -  ^  (uB  ♦  g)  (28) 

where  B  *  is  the  region  in  the  x,y—; plane  corresponding  to  Ba  in 

a  0 

the  X,Y— plane. 

The  boundary  condition  (13)  becomes 

x,y  outside  Ba*,  ivi|r  +  ^  =  0  (29) 

The  condition  of  finite  pressure  at  the  trailing  edge  is  now, 

x  =  »r(y),  ivi  +  ~  is  finite  (30) 

ox 

Finally,  the  condition  of  no  energy  reflection  at  infinity  is  written 
in  the  form 

z-^00,  \)f  S  f(x,y,z)e_iKr  (31) 

where  r2  =  x2  +  y2  +  z^  and  where  f  tends  to  zero  as  z  tends  to 
infinity.2 

^This  ensures  that  0  ~  fei(ajt-Kr) 

as  z  tends  to  infinity 
and  therewith  that  waves  are  traveling  away  from  the  source  of  the 
disturbance . 
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The  boundary  condition  (29)  is  made  more  specific  in  the  following 
manner.  From  condition  (29), 


1r(x,y,0)  =  c(y)e-iVx  (32) 

For  any  line  y  =  Constant  which  does  not  pass  through  Ra*,  it  can  be 
concluded  from  the  condition  of  undisturbed  flow  at  x  =  —00 
that  c(y)  =  0.  The  same  can  be  said  for  that  portion  of  any 
line  y  =  Constant  which  is  situated  in  front  of  the  leading  edge 
of  the  airfoil  region.  The  situation  is  different  for  portions  of 
lines  y  =  Constant  which  do  pass  through  1^*  and  which  are  to  the 

rear  of  the  trailing  edge  of  Ra*.  The  region  to  the  rear  of  the 

trailing  edge  and  bounded  by  lines  y  =  Constant  which  are  tangent 
to  Ra*  is  called  the  wake  region  and  is  designated  by  Rw*»  The 
exterior  of  the  region  Ra*  and  Rw*  is  called  the  remaining  region 
and  is  designated  by  R^.*.  Then  c(y)  =  0  in  Rj,*  and,  in  general, 
c(y)  ±  0  in  Rw*. 

Equation  (29)  is  thus  seen  to  be  equivalent  to  the  following  two 
equations: 


x,y  inside  Rr*,  t  =  0 


(33a) 


x,y  inside  Rw*j  t  =  c(y)e  ivx 


(33h) 


There  is  now  introduced  a  function  A  defined  by 


A(y)  =  2^[»p(y), y,6] 


(3*0 


Equation  (33^)  can  then  be  written  in  the  form 


=  -^A(y)eivM 
dx  2 


x,y  inside  Rw*.> 


(35) 
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In  view  of  equations  (33a)j  there  may  also  be  written 

A(y)  =  2  [**  dx  (36) 

where  xL(y)  is  the  coordinate  of  the  leading  edge  of  the  airfoil 
region  Ra*. 

SUMMARY  OF  THE  EESTJLTAHT  BOUTTOARY-YALUE  PROBLEM 


(31a) 
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Various  quantities  occurring  in  these  equations  are  defined  as 
follows : 


X(x,y)  .  2  (38) 

ox 

The  parameters  k,  M->  and-  V  are  defined  "by  equations  (24),  (22), 

and  (27),  respectively.  The  region  Ra*  follows  from  the  airfoil 

region  Ra  hy  multiplication  with  a  scale  factor  l/h  in  the 

x— direction  and  "by  multiplication  with  a  scale  factor  -  M2/b  in 
the  y-direction.  The  region  1^*  is  the  strip  extending  from  the 

trailing  edge  x  =  xqi  to  x  =  <*>,  and  the  region  Rp*  is  the 

remainder  of  the  x,y— plane. 

The  solution  of  the  boundary— value  problem  is  to  be  used  to 
calculate  the  pressure— change  amplitude  p&  at  the  lifting  surface 

in  accordance  with  the  relation 


Pa  =  -,^2_  ei|ax  (ivt  +  (39) 

b  \ 

which  follows  from  equation  (26) . 

The  solution  of  the  problem  as  summarized  will  be  approached 
through  its  reduction  to  an  integral  equation  for  the  quantity  A  as 
defined  by  equation  ( 38) . 
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AN  INTEGRAL  REPRESENTATION  FOR  TEE  YALTJES  OF  bV/bz 


In  this  section  it  is  proposed  to  derive  a  formula  for  the  values 
of  bty/bz,  for  the  purpose  of  Betting  up  the  basic  integral  equation 
of  the  problem  under  consideration.  To  begin,  results  are  taken  which 
in  essence  are  known  and  then  are  transformed  in  a  way  designed  to 
facilitate  the  subsequent  transition  from  the  exact  double— integral 
equation  of  the  problem  to  the  approximate  integral  equation  containing 
single  integrals  only. 

The  first  formula  is  a  representation  of  the  values  of  dt/cix  in 
the  interior  of  the  half  space'  z  >  0  in  terms  of  the  values  of  bty/bx 
on  the  boundary  z  =  0  of  the  half  space,  as  follows: 

(40) 

In  equation  (kO)  and  in  all  that  follows. 


r2  =  (x  -  |)2  +  (y  -  tj)2  +  z2  (kl) 


and  the  quantity  X  is  according  to  equation  (38)  given  by 
S,ij,0)/SS.  It  is  noted  that  this  representation  of  di/dx 
ensures  that  the  conditions  at  infinity  as  expressed  in  equations  (31) 
are  satisfied. 

From  equation  (kO),  it  follows  that 


If  it  is  now  observed  that  the  quantity  r~ I  e— iKr  is  a  solution  of 
the  differential  equation  (23)  for  t,  then  equation  (k2)  may  be 
written  in  the  alternate  form 
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In  equation  (43),  it  is  noted  that  d2/dy2  =  d2/dq2  and  the  term  in 
question  is  integrated  by  parts  with  respect  to  T).  In  addition  to 
this,  by  making  use  of  the  obvious  identity 

X(  tj)  =  X(t,y)  +  [x(t,q)  —  X(£,yj] 


equation  {43)  is  written  as  follows: 


In  the  first  integral  on  the  right  of  equation  (44)  the  integration 
with  respect  to  q  may  be  carried  out  explicitly.  When  the  remaining 
two  integrals  are  absent,  there  is  thereby  obtained  the  appropriate  form 
of  the  integral  relation  (44)  which  would  follow  if  two-dimensional  flow 
had  been  assumed  from  the  beginning. 

The  following  formula  expresses  the  integral  in  question  in  terms 
of  a  Hankel  function 


— i k\/( x—  | ) 2+  ( y— q ) 2+ z2  (2\  f  [.  T2  2I  ,,  . 

3  dq  =  i«H 0'2)  \_k\JU  -  i)  +  zj  (^5) 


l-oo  V(x  -  |)2  +  (y  -  t])2  +  z2 


The  next  step  is  to  reduce  the  second  integral  on  the  right  of 
equation  (44)  to  an  integral  involving  6x/6q.  To  this  end,  there  is 
written 


[x( !,q)  -  X(l,y)]  9  dq 


(46) 
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The  two  integrals  on  the  right  of  equation  (46)  are  integrated  hy  parts 
and  the  constants  of  integration  are  chosen  in  such  a  way  that  the 
integrated  portions  vanish.  After  some  elementary  transformations, 
this  leads  to  the  following  formula: 


[1(1,11)  -  iU,y)] 


e— 1k\/(x— |)2+(y— T])2+z^ 
f(x  -  i)2  +  (y  -  tj) 2  +  zL 


dTj  = 


Sx  |y  -  tiI 
drj  y  “  1 


dt|  (4 7) 


Equations  (45)  and  (47)  are  introduced  into  equation  (44)  and  the 
following  relation  is  obtained: 


The  final  step  now  consists  in  integrating  equation  (48)  with 
respect  to  x  between  the  limits  —  <»  and  x.  The  condition  of 
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undisturbed,  flow  far  in  front  of  the  lifting  surface  makes  (cM/dz)x_  _ 
There  is  then  obtained  the  relation 


_1_ 

4* 


dX, 

8r] 


r 


^  ly  ~  71 

1  7-1 


s_ 

8x 


0— iK  \J(x—£ )  2+z^+£  2 

\l(x  -  i)2  +  z2  +  £2 


d£  + . 


|y  -  h! 
y  - 


e-itc\/(xt-g)2+z2+^2 
V(x‘  -  |)2  +  z2  +  t2 


dx*  + 


n* 


u 


— 00 


d_ 

q— i  k\Z(x»— 6)  2+z2+(y-q) 2 

\/(xf  -  i)2  +  z2  +  (y  -  q)2 

dx* 


r  d|  dq 


(49) 


Equation  (49)  expresses  the  values  of  cM/dz  for  z  >  0  in  terms 
of  the  values  of  A(|,q)  =  2dt({;,q,0)  /8|  and  is  thus  the  result,  the 
attainment  of  which  was  the  aim  in  the  present  section. 


THE  INTEGRAL  EQUATION  OF  THE  OSCILLATING  LISTING  SURFACE 


The  integral  equation  of  the  problem  is  obtained  by  substituting 
the  information  contained  in  equations  (28),  (32),  and  (35)  in 
equation  (49)  and  by  then  letting  z  tend  to  zero.  In  so  doing  the 
integrals  which  contain  infinities  in  the  integrand  are  to  be 
interpreted  as  Cauchy  principal  values.  This  latter  step  has  been 
discussed  in  detail  for  the  case  of  incompressible  flow  in  reference  1. 
No  additional  difficulties  in  this  regard  are  encountered  for  the 
present  problem  of  compressible  subsonic  flow  and  the  explicit 
Justification  for  this  step  is  therefore  omitted  in  the  present  report. 
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Equation  (50)  holds  when  x  and  y  are  inside  R&*  and  is  to  "be 

solved  for  X,  in  terms  of  A  and  g,  where  A  and  g  are  defined 
by  equations  (36)  and  (37),  respectively. 

For  the  case  of  two-dimensional  flow*  Sl/ch)  =  0  and  the  last  two 
integrals  in  equation  (50)  are  absent.  The  remainder  of  the  present 
work  hag  as  its  object  the  derivation  of  a  procedure  to  take  account  of 
these  last  two  integrals  in  an  approximate  manner  which  permits  the 
calculation  of  the  effect  of  three-dimensionality  of  the  flow  in  a 
simpler  way  than  by  actually  solving  the  complete  equation  (50)*  I*1 

the  derivation  of  this  procedure  the  integral  over  the  airfoil 
region  Ra*  and  the  integral  over  the  wake  region  Ry*  are  treated 

separately. 

REDUCTION  OF  THE  DOUBLE  INTEGRAL  OVER  THE 
AIRFOIL  REGION  Ra* 

Write 
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where 


When  dealing  with  lifting  surfaces  of  sufficiently  elongated 
form,  that  is,  with  surfaces  of  sufficiently  high  aspect  ratio, 
j y  —  rj |  »  |x  —  over  the  major  portion  of  the  surface.  When  this 

latter  inequality  holds,  the  terms  In  ' can,  as  will  he  seen,  he 
written  in  such  a  way  that  a  relatively  simple  dominant  term  can  he 
separated  in  each  of  them.  These  dominant  terms  will  he  used  for  the 
approximation  to  he  developed.  There  remains  then  the  question 
concerning  the  validity  of  the  approximation  over  that  portion  of  R&* 

where  the  inequality  |x  —  £|  «  |  y  —  rj |  does  not  hold.  This  question 
is  answered  as  follows.  It  is  assumed  that  Sk/dq  varies  sufficiently 
slowly  with  T)  so  that  in  this  portion  of  the  region  Ra*  it  is 

effectively  constant.  If  this  is  the  case,  all  that  is  necessary  is  to 
take  account  of  the  fact  that  both  the  kernel  G  and  the  approximation 
to  G  to  he  obtained  are  odd  functions  of  y  —  ij,  80  that  in  both 
cases  the  contribution  to  the  value  of  the  integral  coming  from  this 
portion  of  Ra*  can  be  neglected. 

The  aforementioned  argument  is  also  implicit  in  the' earlier 
derivations  for  the  corresponding  problem  for  incompressible  flow 
(references  1  and  2).  There  appears  to  be  no  reason  to  believe  that 
this  particular  argument  should  be  less  applicable  to  the  problem  of 
subsonic' compressible  flow  than  to  the  problem  of  incompressible  flow. 
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In  equation  (58),  note  that,  when  K  is  not  too  small,  and 
when  |y  -  tj)  »  |x  -  l| , 
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-1k\/(x")2+£2 
V^U")2  +  z2 


(*-« 


(60) 


Neglect  the  term  in  equation  (60)  compared  with  the  term  in 
equation  (59) >  then. 


P-Ml 


Io  Z'** 


^b0<2)(k|?|)  d£ 


(61) 


Note  that  this  approximation  ceases  to  he  correct  as  k  becomes 
smaller  and  smaller.  However,  when  this  is  the  case  the  contribution 
to  the  total  coming  from  I 2  becomes  negligible  because  of  the 

p 

factor  k  in  front  of  the  integral  in  I2. 

Finally,  write  for  I3 


I 


3 


CiT) 


iit 

2 


+ 


Px- 1 


0 


e-iK\/(x,,)2+(y-Tl)2 
V(x  n)2  +  (y  -  tj)2 


(62) 


Now,  when  rcjy  —  T]  is  sufficiently  large. 


2k 
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ft  x  1  -±K\IU"  )2+(y_n) 2  ■  (l  _  6)e-i"|y--i| 

e,.  == - -  dx*  *  =  Qm  1 - - - 5 -  + 

^  +  (v  —  n)2  l.v  —  nl 


/(x*1)2  +  (7  -  T])2 


|y  - 


Rz_-i)K  e-iK[y-nl  (ft) 

Liy  - 


Neglecting  the  terms  in  equation  (64)  compared  with  the  terms  in 
equation  (63), 


i3  a  |;[t  Sb(2)(«|y  "  nl)]  -  f  'yrf- %<2)(*|y  -  n|)  (65) 


Comparison  of  expressions  (57),  (6l),  and  (65)  shows  that  the 
contribution  of  Iq  to  the  total  may  also  he  neglected.  Introducing 

then  expressions  (6l)  and  (65)  into  equation  (53) j  the  following 
approximation  is  obtained 


Equation  (66)  contains  the  fundamental  simplification  of  the  kernel  G- 
in  the  region  Ra*. 
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Noting  that  the  factor  of  dk/ciq  in  equation  (66)  does  not 
depend  on  g,  there  is  written  further 


^(tj) 


dl  =  J- 


^(n) 


fr>*5**&>*5*  (67) 


Take  for  '  (X.)^  its  vai-ue  zero  immediately  in  front  of  the  leading 
edge  and  for  (l)^  its  value  — ivA  which  follows  from  equations  (35) 
and  (38)*  Therewith,  and  with  equations  (36),  equation  (67)  becomes 


\xt(t]) 


^  as  ,  e— iV2T(i)  a_(A -iv*r) 

oq  dq v 


By  introducing  equation  (68)  into  equation  (66),  the  following  equation 
is  obtained: 


&  G  d£  dq  =  1  e-iVST  l_(Ae iv^)  iLjlJlL  ,2 
cSq  J  dq  y  —  q 


■|y-T)| . 


Ho(2)((tKl)d^ 


Ei ^  (K|y  -  n|)|  iTi 


Equation  (69)  represents  the  final  result  of  the  present  section. 


REDUCTION  OF  THE  DOUBLE  INTEGRAL  OVER 


THE  WAKE  REGION  R^* 


In  reducing  the  double  integral  over  the  airfoil  region  Ra*  to 
a  simpler  approximate  form,  use  has  been  made  of  the  fact  that  the  span 
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of  B&*  is  appreciably  larger  than  the  chord  of  Ra*.  Evidently  this 
is  not  the  situation  for  the  wake  region  Rw*  and  thus  additional 
considerations  are  necessary  for  the  integral  extended  over  R^ 

Proceed  as  follows.  Write  the  last  integral  on  the  right  of 
equation  (50)  in  the  form 


The  second  of  the  two  inner  integrals  on  the  right  of  equation  (70) 
may  he  treated  exactly  like  the  integral  over  Ba*.  The  first  of  the 

integrals  has  the  property,  as  will  be  seen,  that  it  is  dependent 
on  x  in  a  simple  explicit  form. 

Taking  the  second  integral  first,  there  results,  in  analogy  to 
equation  (5 6), 


dq  ~ 


dq  y  -  tj 

u 


T?  eo(2)(kKI  ) 


+ 


1  Pxt(t)) 

k  ^  H!(2)(K|y  -  ri|)  e-ivl  dl  dq  (7l) 

_  Jx 

The  remaining  integral  with  respect  to  I  in  expression  (71)  will  be 
introduced  in  evaluated  form  in  the  final  collected  form  of  the  results. 


MCA  TN  1953 


27 


Next  the  remaining'  integral  in  equation  (70)  is  transformed  as 
follows 


dtj 


0  00 


e  1V^X+a^G(-c,  y  -  rj)  dcr 


JO 


dn 


(72) 


It  may  he  noted  that  the  function  F^  reduces  to  the  function  F 

first  introduced  by  Cicala  when  M  =  0,  that  is,  when  the  fluid  is 
incompressible  (references  4  and  5).  Combination  of  equations  (72) 
and  (52)  results  in  the  following  form  of  FM, 


%  = 


e  iV0G(— 0,  y  —  tjj  k)  dcr 


|y  ~  Til 
y  -  t) 


-iva  J 


o 

n_<j 

—CO  \U 


n- 


\™\  1 


a2  +  ^2  2  ~  f 


+  Ik  1  d£  + 


1^1  —lK\jr^+^ 

- 


d£  dT  + 


n-°  -iK\]r2+(y-^r])2 

|y  -  T)|e  w  |y 


T  +  (y  -  t])‘ 


n/t2  +  (y  -  n)2 


+  i  K 


drV  da  (73) 
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Now  equations  (71)  and.  (72)  are  combined  in- accordance  with 
equation  (70)  in  order  to  obtain  the  following  approximation  to  the 
double  integral  over  IV* : 


— (Ae1VXT)e  iV^G  d|  drj  S 
dq 


— iVx 
e 


|^(AeiVxT)FM(y  -  tj)  dq  + 


=-ivxir  _  g-ivx 

Iv 


— (AeiVXT)  ly  ~  HJ 
d^Ae  '  y  -  n 


X 


In 

2 


Bo(2)(d£|)  <1? 


*  ^  Hl(2)(«ly  -  1l) 


dTj 


(74) 
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THE  APPROXIMATE  INTEGRAL  EQUATION  OF  THE  PROBLEM 


Equations  (7^)  and  (69)  are  substituted  into  equation  (50)  and 
the  following  approximate  integral  equation  is  obtained: 


g  U,j) 


1 

2rt 


P^T 

I  X(i,y)K(x  -  I;  k)  dl  + 


n  °° 


g  A(y)e1V^ 


u 


e  iv^K(x  —  k)  dl  + 


Xqi 


e-lM  J 

1 

a  (.c±v*i)  ly  -  hi 

~ 

u  2 

1_|y”H| 

iit 

J 

dq  "  y  —  T) 

K 

1 

2 

— 00 

T  Ho(%l£l)  d£ 


iit  -r,  (2) 


K  ~  H 


(s|y  -  i|) 


£j(a a1V^)F„d,  (75) 


Equation  (75)  represents  in  preliminary  form  the  result  to 'he 
obtained  in  this  report.  This  result  is  reduced  to  a  somewhat  simpler 
form  as  follows. 

Set 


Ae 


lvl*.o 


(76) 


and  introduce  a  new  dimensionless  spanwise  coordinate  y*  defined  by 


7* 


y 

s'Jl  -  M2 


(77) 


30 


NACA  TN  1953 


If  b  represents  the  semichord  at  midspan  and  sb  the  semispan, 
then  it  follows  from  the  definitions  of  x  and  y  in  equation  (19) 
that  the  coordinate  y*  assumes  values  in  the  interval  (—1,1)  only. 

Further,  set 

xL(y)  =  xL*(y*),  xj(y)  =  VKy*)  (78) 

and  note  that  xL*(0)  =  -1,  ^(0)  =  1.  For  a  rectangular  lifting 
surface,  — aj*  =  xj*  =  1  throughout. 


From 'equat ions  (51)  and  (65),  there  follows  for  the  first  integral 
on  the  right  of  equation  (75) 3 
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The  third  integral  on  the  right  of  equation  (75)  becomes, 
with  k  from  equation  (24), 


P 

.  d  (Voivxp)|y  -  n| 

2 

d'r)  y  -  il 

u 

K 

L  - 

■|y-nl 


£  h^KI)  d£ 


-1 


.(2)(«|y  -  ill) 

dT)  = 

p 

n  . —  |y*-n*l 

^  I, 

dn  |y*  -  4*1 

kM 

dry*  y*  —  tj* 

1 

H  1 

2 

— 00 

f  Ho(2)(|?|)  « 


m  i*  (2)/  kaM  .  A 

_  m2  2  1  WTTP  1  7 


dr)* 


(81) 


Finally,  there  is  written  for  the  fourth  integral  on  the  right  of 
equation  (75) j 


|^(Ae1VXr)FM(y  -  4)  dq  = 


1-1 


d^FM 


sk 


/l  -  M2 


(y*  -  4*) 


dq* 


(82) 
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where  the  function  Fj^  as  defined,  try  equation  (73)  3  ma J  he  written 
in  the  following  alternate  form: 
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Finally,  there  are  introduced  dimensionless  coordinates  x* 
and  I*  as  follows: 


x#  =  ^  ~  fo  +  XL) 
*1  ~  XL 


(85) 


so  that  the  interval  x-^  <  x  <  goes  over  into  the 
interval  -1  <  x*  ^  1.  Further, 


-(Xp  +  xL)  =  Ij* 


(86) 


where  h*  is  the  ratio  of  local  semichord  to  semichord  h  at  midspan. 
Finally,  the  following  quantities  are  introduced: 


K*  =  .!** 

V*  =  h*V  > 
ft*  =  $/h* 


(87) 


Note  that,  for  the  rectangular  plan  form,  t>*  =  1,  :%*  =  0  and,  for 

the  elliptical  plan  form,  h*  =  \Jl  —  y*2,  xm*  =  0. 
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Equation  (84)  then  assumes  the  following  form: 


g*(x*,y*) 


_i_ 
2  « 


">1 

:  -  £*)]  d£*  + 


U 


-1 


gi  n*(y*)e-lv*^* 


Hco 


-iV*|* 


K*n[k*(3*  -  I*)]  d{*  + 


3-iV*(x*+Xm*) 
kies'll  -  M2 


4  —  sM 

_1  W 


ks 


_(y*  -  T]*) 


1  -  Mc 


dij*  (88) 


Comparison  of  equations  (88)  and  (84)  gives  for  the  kernels  K 
and  SM  the  following  expressions: 


N(x) 


in  Ixl 
2  * 


%(2)(ui)  45 


(89) 


Sm(x) 


In 

2 


ksM 

)/l  -  M2 


• 

^-M|x| 

|x| 

X 

Hi(2)(M|x|)  + 

Ho<2)(|5|)  4? 

—  00 

i  7==:  Fm(x)  (90) 

Vl  -  M2 


The  taslf  from  here  on  is  the  following.  Equation  (88)  must  be 
solved  for  X*,  in  terms  of  g*  and  ft*.  This  part  of  the  problem 
is  exactly  as  in  the  two-dimensional  theory.  The  function  ft*  is 
then  to  be  determined  by  an  integrodifferential  equation  which  is 
obtained  by  expressing  ft*  in  terms  of  X*  in  accordance  with  the 
definition  of  ft*.  This  part  of  the  problem  is  similar  to  the  earlier 
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work  on  incompressible  flow  of  references  1,  2,  and  6.  Finally  the 
solution  of  the  integral  equation  for  Q,*  must  be  used  to  obtain 
expressions  for  the  pressure  distribution  at  the  airfoil,  as  affected 
by  the  three-dimensionality  of  the  flow  about  a  wing  of  finite  span. 

The  results,  as  expressed  by  equations  (88),  (89),  (90) s  sxlcL  (83) 
include  the  special  case  of  two-dimensional  flow  for  which  dft* / dry*  = 
and  the  special  case  of  incompressible  flow  for  which  M  =  0.  They 
also  include  essentially  known  results  on  compressible  steady  flow  for 
which  k  =  0. 

Massachusetts  Institute  of  Technology 

Cambridge,  Mass.,  May  18,  1948 
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